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C -Abstract 
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Scintillation and ionisation yields for nuclear recoils in liquid xenon above 10 keVm- (nuclear recoil energy) are deduced 
' I / from data acquired using broadband Am-Be neutron sources. The nuclear recoil data from several exposures to two 
D Sources were compared to detailed simulations. Energy-dependent scintillation and ionisation yields giving acceptable 
^ fits to the data were derived. Efficiency and resolution effects are treated using a light collection Monte Carlo, measured 
c/2 photomultiplier response profiles and hardware trigger studies. A gradual fall in scintillation yield below ~ 40keVnr is 
^ ^ found, together with a rising ionisation yield; both are in agreement with the latest independent measurements. The 
^ 'analysis method is applied to the most recent ZEPLIN-III data, acquired with a significantly upgraded detector and a 
O precision-calibrated Am-Be source, as well as to the earlier data from the first run in 2008. A new method for deriving 

the recoil scintillation yield, which includes sub-threshold SI events, is also presented which confirms the main analysis. 
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ZEPLIN-III [H-Q is a dark matter search instrument for 
the direct detection of weakly interactive massive particles 
(WIMPs) via their elastic scattering from xenon target 
nuclei. It records nuclear recoil events via two response 
channels: scintillation and ionisation. A single array of 
photomultipliers (PMTs) records two light signals: one 
prompt, due to scintillation in the liquid (SI); and the 
other delayed, due to the electroluminescence of ionisation 
pharge drifted into the gas region (S2). To use the two 
signals as energy estimators, and so infer the spectrum of 
a WIMP-nucleon scattering population, energy-dependent 
yields of scintillation and ionisation must be established. 
We adopt the conventional definitions: Ces{E) is the scin- 
'tillation yield for nuclear recoils of energy E (commonly 
denoted as true nuclear recoil energy in keVnr) relative 
to that of electron recoils of 122 keV photoabsorption at 
zero electric field; Qy{E) is the yield of ionisation charge 
leaving the interaction site, in electrons per unit energy 
(independently of any other efficiency). 



For a nuclear recoil with energy deposit E in keVnr the 
reconstructed energy (SI) in keVoo is related to CeS by: 



(SI) = ^C,s{E)E, 



(1) 



where 5'nr is the electric field suppression factor for nuclear 
recoils of energy E, and 5*66 is the corresponding factor 
for electron recoils, both defined as unity at zero applied 
electric field. Similarly, using the reconstructed energy 
of the ionisation signal (S2) (as calibrated using electron 
recoils of 122 keV), Qyis defined by: 



(S2) 



W qohm) 



Qy (E) E, 



(2) 
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where W is the energy required to produce an electron-ion 
pair in liquid xenon at infinite electric field (W-vahie) 
and the ratio g(|E|)/go represents the relative fraction of 
charge collected at a finite field. Models and parameteri- 
sations of g(|E|) are taken from Ref. 

The present work measures both yield parameters by 
fitting simulations to observed scintillation and ionisation 
spectra for neutron calibration with an Am-Be {a,n) 
source. This spectral-matching approach has previously 
been applied by others to argon [ij and to xenon f?!, Isj and 
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Figure 1: The discrimination parameter (S2/S1) for the nuclear recoil 
response to an Am-Be neutron source. The colour scale indicates the 
counts for a 4.9 h exposure in the FSR configuration. The analysis 
thresholds for scintillation (SI) and ionisation (S2) are indicated by 
dotted lines. The centroid of the population and the ±2 a contour 
used for this analysis are indicated by the dashed and dash-dotted 
lines. The 40 and SOkeVcc inelastic populations can also be seen. 



was also used by ZEPLIN-III in analysing its first science 
run [9] . Alternatively, measurements with mono-energetic 
neutrons can derive event-by-event recoil energies from 
the kinematics of the scattering angles. Such beam mea- 
surements rely on small prototype chambers rather than 
larger WIMP detectors operating underground. Several 
of these neutron-tagging measurements exhibiting agrad- 
ual fall below 40 keVnr have been reported recently lOj, |ll | . 



We proceed by describing relevant details of event acqui- 
sition and selection from neutron calibration; the Monte 
Carlo simulation of calibration data; estimates of detec- 
tor resolution and efficiency and their application to sim- 
ulated events; and the process of spectral fitting to obtain 
Cce and Qy, including the resultant confidence intervals. 
We include spectral yield measurements from the first Q 
and second science run configurations of ZEPLIN-III, hav- 
ing re-analysed the earlier data with improved software. 
The results are further consolidated by using an indepen- 
dent new method to analyse events below the scintillation 
threshold. £off is therein derived by attributing the aver- 
age scintillation signals to the nuclear recoil energies re- 
constructed using the ionisation yield Qy. The combined 
effects of the revised efficiency and light yield on WIMP- 
nucleon scattering limits are presented. 

2. Neutron data 

The first science run (FSR) configuration and data sets 
of ZEPLIN-III are described in detail in Ref. ;93. The 
main neutron calibration data set was acquired in May 
2008. For the 2010-11 second science run (SSR) the PMTs 
were replaced by new units, with a 40-fold reduction in 
radioactivity [12i| and an anti-coincidence (veto) detector 



was installed around the experiment ISj, 1141 . Other 
changes include the use of a calibrated neutron source 
positioned centrally above the active volume and a ~ 10% 
lower drift field in the liquid of 3.4kVcm^^. The trigger 
for event acquisition was obtained from a shaped sum of 
PMT signals which, at low energies, was derived from the 
S2 signal. 

In the SSR, the Am-Be source was inserted through 
the experiment shielding via a dedicated delivery pipe, 
coming to rest above the cryostat vessel just ~ 5 cm offset 
from the PMT array centre. The software acceptance 
threshold for SI pulses requires a 3-fold coincidence 
amongst the 31 PMTs of at least one photoelectron in 
each channel. In this analysis, only events representing 
single elastic scatters of neutrons are considered. The 
event selection is achieved using only basic waveform 
and pulse quality cuts. Interactions of photons from 
the source, inelastic scattering and radiative capture on 
xenon and other materials are rejected by event selection 
in S2/S1, as shown in Fig.[T] Further selection criteria, 
e.g. to remove MSSI events (see Ref. ;9J for details), have 
been found to have no significant effect on the results of 
this analysis. The SSR neutron data were acquired with a 
calibrated 20MBq Am-Be source emitting 1321±14n/£[3 
for a total of 10.2 h on three different days between June 
and August 2010. A detailed analysis of the effect of the 
operational parameters showed no significant difference 
between the various datasets. The primary scintillation 
spectrum of single elastic scatters within a 3.5 kg fiducial 
volume in the centre of the detector is shown for both 
FSR and SSR in Fig.m 

The SSR source was used to calibrate that employed in 
the first run using a low-background HPGe detector, re- 
sulting in 5512 ± 358 n/s. In the FSR calibration (5-hour 
exposure) the neutron source was only marginally inserted 
into the shielding to avoid too high an event rate and ir- 
radiated the detector from one side, which changes the 
expected recoil spectrum slightly. 

3. Neutron scattering simulation 

Geant4 15, [l3| Monte Carlo simulations of the entire 
ZEPLIN-III experiment have been implemented to predict 
the detector response to neutrons [2,]. For the simulations 
in this work, Geant4.9.3. and updated Xe(n,n) scattering 
cross sections from the ENDF /B-VII database were 
used. A comprehensive and accurate setup of the exper- 
iment was implemented, taking into account the detector 
configuration and source position for the two data sets. 
The simulations include the liquid and gaseous xenon, the 
PMT array, electrode grids, copper vessels, the source de- 
livery mechanism, the plastic and lead shielding and, for 



t Calibration by National Physical Laboratory, Teddington, Mid- 
dlesex, UK, May 2009 
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Figure 2: The differential spectrum of the primary scintillation of 
single scatter nuclear recoils from the 2008 (top) and 2010 (bottom) 
Am-Be data. Also shown are the results from a detailed Monte Carlo 
simulation (solid lines) considering an energy dependence of C^s as in 
Fig. |3] and the corresponding nuclear recoil energy scale. SI energies 
below 2.0keVcc(FSR) and 2.5keVco(SSR) are not used to obtain the 
fit and are shown as extrapolations only (dotted). Superimposed is 
the detection efficiency (dash-dotted). 



light collection efficiency. To study the detector response 
and assess the evolution of a multitude of operational 
parameters over the complete data-taking period, an 
external ^''Co source was inserted into the shielded vol- 
ume every day. The 122 keV interactions, together with 
Compton-scattering events observed during calibration 
with a ^'^^Cs source, are used to determine radial response 
profiles for all channels for position reconstruction in 
the horizontal plane (both in SI and S2). As a result, 
the average light detection efficiency at the operating 
field is 1.8 ±0.1 and 1.3 ± 0.1 phe/keVco for the FSR 
and SSR, respectively. Applying the response profiles 
for each PMT to Monte Carlo simulations of nuclear 
recoils, one can determine the fraction of events in each 
bin of reconstructed energy meeting the required 3-fold 
coincidence. Alternatively, relaxing the coincidence 
requirement to 2-fold and studying the fraction of 2-fold 
to 3-fold events in each bin and comparing with Monte 
Carlo simulations gives a similar detection efficiency. The 
resulting detection efficiencies are shown in Fig. [51 

The daily calibrations are also used to determine further 
correction factors. For example, the depth-dependence 
of the ionisation signal S2 varies with the mean electron 
lifetime in the liquid, which must be corrected accurately. 
The width of the electroluminescence signal is propor- 
tional to the thickness of the gas phase and hence its 
polar distribution depends on the detector tilt, which has 
been observed to vary slowly due to geological factors. 
Additionally, the electroluminescence yield varies with gas 
pressure. All of these corrections are taken into account. 



the SSR configuration, also the veto detector. The simu- 
lated energy spectrum for neutrons from the Am-Be source 
extends up to llMeV with a mean of 4.1MeV [l3|. The 
low-energy threshold was 100 keV; no significant spectral 
uncertainty is to be expected above this energy for the 
low activity source considered here. That corresponds to 
a maximum recoil energy of 3 keVnr in xenon, which is sig- 
nificantly lower than the analysis threshold, and so this 
uncertainty has limited impact on the results. Previous 
studies showed the Monte Carlo results for low energy nu- 
clear recoils to be relatively insensitive to variations in the 
neutron source spectrum, the precise source location and 
the effect of intervening and surrounding materials 9] . The 
Monte Carlo event selection includes nuclear recoils from 
elastic and inelastic scattering (when prompt 7-rays may 
escape from the xenon undetected); recoils from radiative 
capture (when 7-rays can be delayed significantly) affect 
the spectrum only below 1 keVnr and were not considered. 

4. Detector response 

The overall detection efficiency is dependent on the 
quantum efficiency of the photomultipliers, on hardware 
and software trigger efficiencies and on the geometrical 



The scintillation yield of liquid xenon depends on the 
applied electric field. The recombination of ionisation 
is partly suppressed by the electric field, resulting in a 
reduced scintillation response. The scintillation quenching 
for electron recoils, Soo, has been measured in ZEPLIN-III 
(19I I and found to agree with other published data pol . [2ll | . 
In the FSR configuration, Soo = 0.38; for the slightly 
lower field configuration during the SSR, Soo = 0.39. The 
field-induced quenching factor for nuclear recoils has been 
measured by other experiments and is found to neither 
change significantly with field nor be very dependent on 
the nuclear recoil energy 11, |2^. The value for Snr is 
taken to be 0.92. 

According to Eq.[51 the total number of electrons ex- 
tracted into the gas phase for a given signal S2 can be de- 
termined with the help of the detector response to single 
electrons extracted from the liquid. The signal size of sin- 
gle electrons, SE, has been measured consistently and with 
high efficiency from two types of data: random triggers 
and searches for photoionisation between SI and S2: SE = 
30.6 ± 0.5phe/e^ in the FSR and SE = 11.8 ± 0.4phe/e- in 
the SSR [2^ [23 1 . The electron emission probability at the 
hquid/gas interface is 77 = 0.83 (FSR) and 0.66 (SSR)!!!]. 
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Figure 3: The energy-dependent relative scintillation yield for nu- 
clear recoils Z^eff (solid lines) for the two ZEPLIN-III datasets, in- 
cluding relevant 68% C.L. bands (FSR: green \\, SSR: blue //). Below 
the analysis range (corresponding to ~ 7-9keVnr), the scintillation 
yield is indicated as dashed lines. Also shown are previous pub- 
lished measurements using mono-energetic neutron beams: (>) llClll . 
(■)[I1. (^)[11. (0)111, ('^)[27l, (v);28] and (O) 21^, or obtained 
using a similar Monte Carlo matching procedure (t)i7|1. 

5. Fitting and uncertainties 

5.1. Relative scintillation efficiency 

The nuclear recoil energy spectrum obtained from the 
Monte Carlo simulations (in keVnr) can be converted to 
a spectrum in reconstructed energy (SI) (in keVoo) via 
Eq.[TJ At low energy, the energy resolution applied to the 
Monte Carlo is defined by Poisson fluctuations in the num- 
ber of photoelectrons and PMT gain fluctuations and was 
modelled by a continuous Poisson distribution replacing 
the factorial by Euler's F function. The resulting con- 
volved energy spectrum is multiplied by the overall detec- 
tion efhciency and compared with the experimental Am-Be 
data (see Fig. [5]). To allow optimal freedom and avoid any 
model bias, natural piece-wise cubic splines with continu- 
ous 1st and 2nd derivatives are used to parameterise the 
energy dependence of £cff , a method similar to that used in 
Ref. [7| ■ The global best fit to the data in the energy range 
between 2.0keVoo(FSR), or 2.5 keVoe(SSR), and lOOkeVoc 
is found by using a minimum technique and is shown in 
Fig. [3] The lower analysis thresholds are determined by an 
approximately 50% detection efficiency. In this work, the 
spline points are fixed at 0.5,2.5,6.0,15,50 and 200keVnr 
and are unconstrained in £off- The outcome of the fit has 
proven to depend very little on the position and number 
of spline points (since the best fit functions are almost fea- 
tureless) . The dashed lines in Fig. [2] show the differential 
rate as a function of the reconstructed energy of the Monte 
Carlo data using the derived form of £eff for the two data 
sets. 

5.2. Recoil ionisation yield 

A similar approach to determine Ccs, is also applied 
to the energy dependence of Qy. In general the W- value 
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Figure 4: The differential spectrum of secondary scintillation signals 
from single scatter nuclear recoils from the two datasets (top: FSR; 
bottom: SSR). Also shown are the Monte Carlo data (solid lines) con- 
sidering an energy dependence of Qy as in Fig.[5]and the correspond- 
ing nuclear recoil energy scale. S2 energies below 1.0keVoc(FSR,) or 
f .2keVcc(SSR) and above llkeVeo are not used to obtain the fits 
and are shown only as extrapolations (dotted). Superimposed are 
the detection efficiencies (dash-dotted) as obtained by the primary 
scintillation signal SI (see Fig.O and converted into S2 energy using 
the mean correlation between the two channels. 

is defined as the average energy required to produce an 
electron-ion pair. Using the number of electrons recorded 
in S2 for 122keVco 7-rays from ^^Co, we can confirm 
the W-value for liquid xenon by rearranging Eq.[21 For 
SSR data we obtain W = 16.5 ± 0.8 eV, which agrees very 
well with the reference measurement of W = 15.6 ± 0.3 eV 
published in Ref. [131 ■ ^ rate-dependent photocathode 
charging effect for the FSR PMTs prevented an accurate 
W-value being derived from the FSR data. This effect has 
not been seen in the primary scintillation signal SI, and 
hence is irrelevant for the determination of the energy 
scale for low energy recoils and does not infiuence the 
results reported in Ref. For the following analysis 
of the ionisation yield and the comparison to the SSR 
measurement, the W-value was set to W = 15.6eV for this 
dataset. 

The observed nuclear recoil spectra for Am-Be neutrons 
are converted to number of electrons via Eq. [5] and con- 
volved using Gaussian deviates into a number of electrons 
emitted from the liquid. Applying the same technique as 
for yCeff , the obtained spectra are then matched to the data 
by varying the energy dependence of Qy. Spline points are 
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Figure 5: The ionisation yield Qy for nuclear recoils as derived from 
the FSR (green \\) and SSR (blue //) datasets including relevant 68% 
C.L. bands. Results below the analysis range are indicated by dashed 
lines. Also shown are previous measurements at l.OkV/cm (#) and 
4.0kV/cm (■) from Ref. O, at 2kV/cm (□), (A), 0.3kV/cm (O) 
and O.lkV/cm (O) from Ref. [2ll and spectra obtained using similar 
Monte Carlo matching procedures at 0.73kV/cm (A) 7] and (▼) [29| . 
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Figure 6: The nuclear recoil response to Am-Be neutrons in the SSR 
configuration including events below the previous analysis threshold 
for SI and S2 (dotted lines). The colour scale gives the counts for a 
5.9h exposure. Also indicated are the ±3cr contours to the centroid 
of the neutron recoil population (dashed and dashed-dotted) and 
their extension sub-threshold. 



fixed at 0.5,4.0,10,30,75 and 250keV„r. Fig.H shows the 
differential rate in ZEPLIN-III and the relevant detection 
efficiency for the FSR and SSR data. In contrast to the 
scintillation signal, the overall detection efficiency for the 
ionisation signal is predominantly determined by SI, since 
the analysis of data requires a primary signal. For high 
nuclear recoil energies, both SI and S2 are detected with- 
out any efficiency losses. Below SI signals of ~5keVoo, 
corresponding to S2 signals for nuclear recoils of ~2 keVco, 
the correlation between ionisation and scintillation signal 
becomes unreliable due to the low number of SI photo- 
electrons. Hence, to avoid any dependence on the effective 
scintillation yield, CcS, the ionisation yield is only deter- 
mined for nuclear recoil energies above lOkeVnr and the 
extension to lower energies is indicative only. The result is 
shown in Fig.[5l together with results of a similar approach 
by XENONIO 0,l3, as well as measurements at various 
electric fields and nuclear recoil energies using fixed neu- 
tron energy scattering experiments 111. |21|. 



5.3. Alternative method to determine Cce 

In order to provide a cross-check to the scintillation 
yield, Ccs, a further analysis has been developed omitting 
the dependency on the scintillation detection efficiency 
on the Monte Carlo simulations and relying instead on 
ionisation yield only. This has been achieved by relaxing 
the event selection criteria to include any events with an 
ionisation signal and, hence, lowering the energy thresh- 
old for nuclear recoils and avoiding detection efficiency 
penalties in SI. The scintillation signal for a given nuclear 
recoil energy in this analysis is estimated from the average 
amount of light preceding the ionisation signal, regardless 
of the 3-fold coincidence requirement. This results in 
an effective threshold for the scintillation signal close 



to zero. Assuming a power law fit to the experimental 
data of Fig.O the number of electrons of the ionisation 
signal can be converted to the corresponding nuclear 
recoil energy using Eq.[2] The effective scintillation light 
yield is then determined via Eq.[TJ In Fig.[S] the event 
selection for nuclear recoils in this additional analysis 
is shown, including sub-threshold events. Here, only 
events within 120 mm from the centre of the detector are 
considered to avoid any bias due to lower light collection 
of peripheral events. Due to the lack of the primary 
scintillation signal and, hence, resolution in the z-axis, 
the correction due to the mean electron lifetime in the 
liquid is a depth-averaged value for low energy nuclear 
recoils. This is justified by the reasonably homogenous 
distribution of these recoil events in the xenon volume. 



To avoid bias from PMT thermionic photoelectron emis- 
sion, cross-phase single electron emission and other spuri- 
ous effects producing non-correlated VUV light, an average 
background is determined by integrating all pulses over the 
same time period before primary scintillation signals (SI). 
This background was found to be 3.1 ± O.lphe/event dur- 
ing the Am-Be calibration in the SSR. The results of this 
alternative analysis method are shown in Fig.[7l The er- 
rors arise predominantly from the Poisson fluctuations in 
number of photoelectrons and ionisation electrons, as well 
as from the light collection efficiency and include the effect 
due to the fitting of a power law for Qy. To summarise, 
this method was used to derive the effective scintillation 
yield, Ccs, without relying on the Monte Carlo simulations 
or being limited by detection efficiencies and the results are 
in agreement with the spectrum fitting analysis method. 
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Figure 7: The relative scintillation yield for nuclear recoils C^ft as 
determined by extending the event selection to S2-only events for 
the SSR Am-Be data. Also shown are the results obtained using the 
Monte Carlo fitting method, as in Fig.|3] 



5.4- Error analyses 

The 68% confidence intervals, as shown in Fig. [3] and 
O are determined by models of CcS and Qy giving 
< Xmin + Qi' where = 7.01 for six free parameters 
3l| . Another uncertainty beyond the fitting procedure 



contributing to the results of the scintillation and ioni- 
sation yields arise from the detector dead-time and the 
Am-Be neutron source strengths. Both of these param- 
eters scale directly the simulated recoil spectra. The 
neutron rates from the two sources are known within 6.5% 
(FSR) and 1.1% (SSR) and the dead-time of the system 
produces an uncertainty of ~ 2% in the live times of the 
neutron exposures for the two runs. These uncertainties 
propagate into <5% change of CcS and Qy over the 
analysed energy range. 

The results are also affected by systematics of the 
detector and our knowledge of the properties of liquid 
xenon. For example, the photoelectron yield of the 
chamber, Ly, varies with the position of the interaction 
in the detector and is determined by ^^Co calibration for 
122 keV 7-rays. The light collection decreases by ~ 10% 
of the volume-averaged value for interactions more than 
100 mm away from the centre of the detector. However, 
including more peripheral interactions in this analysis 
(resulting in a lower average yield, but more events) did 
not change the results for the effective scintillation yield 
£cffby more than a few percent. These are the dominant 
systematic uncertainties in this analysis. Their combined 
magnitude does not contribute beyond the statistical 
error estimation from the fitting procedure as shown in 
Fig. [31 Other systematic errors, also affecting the beam 
measurements, are discussed critically in Ref. [32| . 



of the fitting procedure. Here, the detection efficiency as 
shown in Fig|31 is afflicted by the aforementioned correla- 
tion between SI and S2 below -lOkeVnr and the signif- 
icant relative variance of photoelectron numbers at these 
energies. Further, the emission efRciency of electrons from 
the liquid into the gas phase is strongly dependent on the 
electric field in the liquid and thus, in the case of ZEPLIN- 
III on the precise thickness of the liquid volume. A con- 
servative systematic error of 5% on 77 directly translates 
to the ionisation yield Qy via the W-value. In the FSR, 
the uncertainty on the used measurement of the W-value 
(2%) from Ref. fsO] has been added to the systematics. For 
Qy, the combined systematic errors exceed the statistical 
errors from the fitting procedure and hence are included 
in the error band in Fig.O 

6. Conclusions 

We have presented the energy dependence of the 
scintillation and ionisation yields of nuclear recoils 
above 10keV„r as seen in the ZEPLIN-III experiment 
using two data sets. The neutron calibration data 
obtained during the first science run (FSR) in 2008 were 
reanalysed using the same pulse finding and position 
reconstruction algorithms as used on the 2010 data. 
The reanalysis of the FSR data revealed a very small, 
but critical population of events whose analysis relied 
on the dual-range data acquisition cross-over point. 
This explains the new £eff curve. The WIMP-search 
results are only marginally affected with now 5 events 
in the (wider) 7-36keVnr signal acceptance region. This 
corresponds to an electron recoil leakage of 1 : 7800 in the 
worst case. The 90% C.L. hmit on the WIMP-nucleon 
scalar cross-section improves for WIMP masses below 
30GeV/c^ and increases for higher masses relative to 
that published in Ref. [9] . For example at the bottom of 
the sensitivity curve (~ 50 GeV/c^) the limit is 23% higher. 

Using a complete new set of PMTs in the second science 
run configuration (with a 40-fold reduction in radiological 
background and very different optical and electrical per- 
formance), a slightly lower external electric field and a 
different neutron source strength and location, consistent 
results were obtained for the scintillation yield. System- 
atic errors resulting from the detector configuration and 
the applied fitting method are understood and show only 
an overall shift within a few percent, and result in little 
change of the spectrum. The results below lOkeVnr are af- 
fected by the overall detection efficiency, the average light 
yield and light collection efficiency of the ZEPLIN-III de- 
tector and should be treated as extrapolations only. 
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